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ABSTRACT 

We analyze the X-ray spectra and their timing properties of the compact X- 
ray binary 4U 1820-30. We establish spectral transitions in this source seen with 
BeppoSAX and the Rossi X-ray Timing Explorer ( RXTE ). During the RXTE 
observations (1996 - 2009), the source were approximately ~ 75% of its time in 
the soft state making the lower banana and upper banana transitions combined 
with long-term low-high state transitions. We reveal that all of the X-ray spectra 
of 4U 1820-30 are fit by a composition of a thermal (blackbody) component, a 
Comptonization component (COMPTB) and a Gaussian- line component. Thus 
using this spectral analysis we find that the photon power-law index T of the 
Comptonization component is almost unchangeable (r ~ 2) while the electron 
temperature kT e changes from 2.9 to 21 keV during these spectral events. We also 
establish that for these spectral events the normalization of COMPTB component 
(which is proportional to mass accretion rate M) increases by factor 8 when kT e 
decreases from 21 keV to 2.9 keV. Before this index stability effect was also found 
analyzing X-ray data for ^-source GX 340+0 and for atolls, 4U 1728-34, GX 3+1. 

Thus, we can suggest that this spectral stability property is a spectral signature 
of an accreting neutron star source. On the other hand in a black hole binary T 
monotonically increases with M and ultimately its value saturates at large M. 
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(4U 1820-30) :radiation mechanisms: nonthermal — physical data and processes 
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1. Introduction 

Accreting neutron stars (NS) can be observationally classified using a color-color dia- 
gram (CCD) into two distinct categories, atoll and Z sources, based on their different CCD 
forms when the source undergoes the spectral and luminosity changes. Along with this 
phenomenological difference between atolls and Z-sources in terms of the CCD there are 
important X-ray spectral and timing characteristics, which are essentially different for these 
types of NS low mass X-ray binaries (LMXB). The main observational difference between 
these types is the specific range of luminosity changes. The atolls are observed when their 
luminosity changes from 0.01 up to 0.5 of the Eddington limit L Edd (see Christian & Swank 
1997; Ford et al. 2000) while Z-sources are seen when the resulting luminosity is near Ed- 
dington regime (e.g. Seifina, Titarchuk & Frontera 2013). In this Paper we present our 
analysis of peculiar atoll 4U 1820-30, which is a bright atoll source in terms of its luminosity 
and, at the same time, it is a typical atoll source in terms of timing evolution. On the other 
hand, during bright phase 4U 1820-30 is as bright as a subclass of persistent bright atolls 
(GX 13+1, GX 9+1, GX 9+9 and GX 3+1), but 4U 1820-30 shows larger range of luminosity 
and demonstrates all states in terms of CCD (from the island to banana states), whereas 
these aforementioned bright atolls are only seen in the banana state (e.g. Hasinger & van 
der Klis 1989). 

Furthermore, 4U 1820-30 has a maximal luminosity about 0.5 of L Edd and thus it adjoins 
to Z source luminosity range. Also 4U 1820-30 shows Type-I X-ray bursts and characteristic 
timing features as a typical atoll demonstrating an evolution of band-limited noise (BLN), 
very low frequency noise (VLFN) components and low frequency quasi-periodic oscillations 
(LFQPOs) (in 20-40 Hz range) when it evolves from the banana to island states. However, 
opposite to the other atolls, 4U 1820-30 exhibits LFQPOs with frequencies near 6 Hz during 
the banana state (Wijnands, van der Klis, & Rijkhorst 1999; also see Sect. 6.2 in this Paper), 
which are usually seen in Z-sources during the Normal branch. Thus, 4U 1820-30 com- 
bines properties of ordinary and bright atolls and Z-sources in terms of timing and spectral 
evolution, luminosity and detection of 6 Hz QPO. 

4U 1820-30 is a LMXB observed at 0.66” from the center of the NGC 6624 cluster. 
Grindlay et ah (1976) were the first who identified this source as a Type-I X-ray burst. 
Kuulkers et al. (2003) estimated the distance d = 7.6 ±0.4 kpc to 4U 1820-30 assuming that 
the peak luminosity equals to L Edd for He burst atmosphere. Vacca, Lewin & van Paradijs 
(1986) estimate the distance to 4U 1820-30 as 6. 4+0. 6 kpc using the analysis of UV diagrams 
for NGC 6624. Rappaport et al. (1987) find that the binary system comprises of a He white 
dwarf of mass of 0.06+0.08 M 0 and NS, [with a mass later evaluated by Shaposhnikov & 
Titarchuk (2004) as ~1.3 M 0 ], orbiting at the period of 11.4 minutes (Stella et al. 1987). 
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Hasinger & van der Klis (1989) classify 4U 1820-30 as an atoll source. Priedhorsky & Terrell 
(1984), Simon (2003) and Wen et al. (2006) find that the flux varies between the soft and 
the hard states ( banana and island ones respectively) are quasi-periodic with period at ~170 
d and these flux variations have been suggested related to tidal effects of a remote third 
star [Chou & Grindlay (2001) and Zdziarski et ah (2007)]. RXTE observations revealed 
4U 1820-30 as a prominent source of kilohertz quasi-periodic oscillations (kHz QPO) (Smale 
1997). X-ray bursts are only observed at low fluxes [e.g. Clark et al. (1977)]. Furthermore, 
Cornelisse et al. (2003) and Zhang et al. (1998) find that the observed kHz QPOs corre- 
late with the flux which probably suggests that these variations are because of a luminosity 
change caused by change of mass accretion rate. 

Strohmayer & Bildsten (2004) established that the short (~10 - 15 s) Type-I outbursts 
are because of the unstable thermonuclear burning of mixture of hydrogen and helium at 
the NS atmosphere bottom. SAS-3 observations showed strong evidence that X-ray bursts 
can only occur in its low-intensity state (Clark et al. 1977). All bursts observed from 
4U 1820-30 indicate that a photosphere expends with an increase of a photospheric radius 
by a factor of 20. Such an expansion leads to strong softening of the resulting spectrum 
(see e.g. Strohmayer & Bildsten 2004). Moreover, a several hour long ’’superbursts” was 
observed from 4U 1820-30 on September 9, 1999. It is now understood that superbursts can 
be caused by the burning in the carbon ashes produced by Type-I bursts (Strohmayer & 
Brown 2002). 

Einstein, EXOSAT, Ginga, ASCA, and BeppoSAX also observed 4U 1820-30. Many 
spectral models have been applied to fit the observed X-ray spectra. For example, models 
which are a sum of a blackbody (BB) with thermal bremsstrahlung or a blackbody with a 
power law combined with exponential cutoff (CPL). A more detailed model were developed 
using a Comptonization spectrum by Sunyaev & Titarchuk (1980) (see CompST model in 
XSPEC) combined with a blackbody. Note, White et al. (1986) and Christian & Swank 
1997) show that the models based on thermal bremsstrahlung mechanism are unphysical. 

The emission measures found using this model have to be of order ~ 10 60 cm -3 and they 
are too large for the plasma cloud near NS which radius is only of order of 10 6 cm. 

Therefore the CPL and CompST components combined with a blackbody were applied 
to fit X-ray spectra of 4U 1820-30 [see e.g. Bloser et al. (2000)]. The BB temperature kT BB , 
the photon index T and exponential cutoff energy Eq are the CPL+BB model parameters. 
Bloser et al. (2000) also included photoelectric absorption at low energies using the cross 
section of Morrison & McCamman (1983). White et al. 1986 and Hirano et al. (1987) show 
that a Gaussian at ~6.7 keV is often needed to take into account a blend of K a iron lines. 
Parsignault & Grindlay (1978) applying a power-law fit to the 4U 1820-30 ANS data (ANS 
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is an abbreviation of the Astronomische Nederlandse Satelliet) found X-ray spectral changes 
due intensity variations. They obtain that the photon index T changes from 2 at high count 
rates to 1.4 when count rate is low. In other words, they find that spectrum becomes harder 
when luminosity decreases. Stella, White, & Priedhorsky (1987) used the CPL+BB model 
to fit the data of the EXOSAT ME instrument in the energy range from 1 to 30 keV. These 
particular EXOSAT data were obtained during 1984 - 1985. The source was found at a 
high luminosity state (6.0xl0 37 erg s~ 3 ) and a low luminosity state (2.0xl0 3 ' erg s _1 ). The 
best-fit parameters of these EXOSAT spectra T, Eq and kTsB are at 1.7, 12 keV and 2 keV 
respectively in the high state while P ~2.5, Ec >30 keV, and kTsB= 2-3 keV in the low 
state. 

Smale et al. (1994) analyzed the ASCA/GIS data for 4U 1820-30 which was observed in 
the low state in 1993. They fit the 0.6-11 keV spectrum using the CompST + BB model. The 
best-fit CompST parameters were around 3.6, 13.5 and 0.76 keV for the plasma temperature, 
optical and a blackbody temperature respectively. Christian & Swank (1997) reported on 
the Einstein (SSS + MPC) 1978 observation in the 0.5 - 20 keV energy range. They found 
the source was in the high state characterized by the luminosity of 5.5 xlO 37 erg s _1 and the 
best-fit parameters of the CompST+BB model were very similar to those obtained using the 
ASCA data. 

Piraino et al. (1999) and Kaaret et al. (1999) analyzed the first observations of 4U 1820- 
30 extended above ~30 keV. For this analysis they used the observations implemented by 
the NFIs (0.1 - 200 keV) of BeppoSAX in 1998. The best-fit of the observed spectrum in the 
0.3 - 40 keV energy range give kT B B= 0-47 keV, T=0.55, and E c = 4.5 keV for the CPL+BB 
model , whereas the best-fit parameters of the CompST+BB model are &Trb= 0.46 - 0.66 
keV, kT e = 2.83 keV and r=13.7. It is worth noting that the instruments with a response 
below 1 keV provide low values of kTsB ■ Although, the range of luminosities and the best-fit 
parameters inferred using BeppoSAX and the CompST + BB model gave very similar values 
with respect to that obtained using the other instruments (see above). 

The BeppoSAX Phoswich Detection System (PDS) could not detect the emission from 
4U 1820-30 above 40 keV. Note, BATSE [see Bloser et al. 1996] also failed to find this 
source in the 20 - 100 keV range during the first four years of the CGRO. For the first 
time a high energy tail above 50 keV has been found by INTEGRAL in the hard state of 
4U 1820-30 (Tarana et al. 2006), which put this source in the list of X-ray bursters which 
exhibit high-energy emission. 

It is interesting to note that also other atolls can be described by similar spectral models. 
Lin et al. (2007), hereafter LRH07, pointed out that at the low-L^ end of the soft-state track 
a weak Comptonization component is needed. LRH07 modified the BPL and COMPTT 
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models applying their modifications to atolls Aql X-l and 4U 1608-52. They were interested 
to find out how much energy is directly visible as a pure thermal radiation and thus one can 
obtain the remaining fraction for the Comptonized radiation /. In this sense this approach 
is similar to that using the COMPTB model (see below §3) . In this way, to account for 
specific spectra LRH07 describe the hard state by a BB+BPL model and the soft state by 
means of a three component model, MCD+BB+CBPL, where CBPL is a broken power law 
with the high energy cutoff taking into account the Comptonization effect. 

In this Paper we show a thorough X-ray spectral-timing analysis of the data for 4U 1820- 
30 using the BeppoSAX and RXTEj PCA/HEXTE available observations which were made 
during 1998 - 1999 and 1996 - 2009 years respectively. Unlike the past spectral analysis, we 
adopt an unified model, capable describing the spectra observed during both the soft and 
hard states. The full list of observations used in our data analysis is present in §2 and Tables 
1 and 2 while we describe, in detail our spectral model and spectral analysis using this model 
in §3. We interpret X-ray spectral-timing evolution when the source undergoes the spectral 
state transition in §§4—6. We explain our results in detail and come to the final conclusions 
in §§7—8. 


2. Data Selection 

We obtain broad band energy spectra of the source using data from three BeppoSAX 
Narrow Field Instruments (NFIs), namely the Low Energy Concentrator Spectrometer (LECS) 
with the 0.3—4 keV energy band, the Medium Energy Concentrator Spectrometer (MECS) 
with the 1.8—10 keV band and the Phoswich Detection System (PDS) with the 15-200 keV 
band [see Parmar et al. (1997); Boclla et al. (1997); Frontera et al. (1997) respectively]. 

We use the SAXDAS data analysis package for data processing. We renormalized the 
LECS data based on the MECS data. We treat relative normalizations of the NFIs as free 
parameters when we proceed with model fitting, but we fix the MECS normalization at 
1. Each of these normalizations is controlled if they are in a standard range for a given 
instrument 1 . Furthermore, we rebinned the spectra to obtain significant data points. The 
LECS spectra are rebinned using a binning factor which varies with energy (Sect. 3.1.6 of 
Cookbook for the BeppoSAX NFI spectral analysis) implementing rebinning template files 
in GRPPHA of XSPEC 2 . The PDS spectra are rebinned with a linear binning factor 2, 


1 http://hcasarc. nasa.gov/docs/sax/abc/saxabc/saxabc.html 

2 http://hcasarc. gsfc.nasa.gov/FTP/sax/cal/responses/grouping 
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namely we group two bins together leading to the bin width of 1 keV. For all of these spectra 
we use a systematic error of 1%. The BeppoSAX observations implemented in our analysis 
are shown in Table 1. 

We also use publicly available the RATE data sets (Bradt et ah 1993) which were 
obtained from April 1997 to March 2009. In total, they include 92 observations taken at 
different states of the source. The LHEASOFT/FTOOLS 5.3 software package were applied 
to process the data. Also for our spectral analysis we apply PCA Standard 2 mode data, 
collected in the 3-20 keV energy range and the most recent release of PCA response 
calibration (ftool pcarmf vll.l). We use the standard dead time correction to the data. 

A background corrected in off-source observations is subtracted from the data. We use 
only data from 20 to 150 keV energy in order to avoid the problems related to the HEXTE 
response and background determination. We apply the GSFC public archive to analyze all 
available data sets (see http://heasarc.gsfc.nasa.gov). We present a full list of observations 
covering the source evolution during different spectral state events in Table 2. 

We implement an analysis of thirteen years RXTE observations of 4U 1820-30 for 7 
intervals (see blue rectangles in Figure 1 in Titarchuk, Seihna & Frontera (2013), hereafter 
TSF13). We fitted the RXTE energy spectra using XSPEC astrophysical fitting software. 
For our data analysis we have also applied the public available 4U 1820-30 data in the energy 
range from 2 to 10 keV from the All-Sky Monitor (ASM/ RXTE) for all observation scans. 

According to ASM monitoring system 4U 1820-30 shows long-term variations with pos- 
sible period ~176 days of the 2-10 keV flux [see Fig. 1 in TSF13 and Priedhorsky & Terrell 
(1984); Simon (2003) and Wen et al. (2006)]. The count rate changes in the interval of 
5—35 counts s” 1 throughout each cycle (see in Fig. 1 in TSF13). Our RXTE spectral stud- 
ies are directed to investigate: i) the continuum spectrum, in particular, the hard X-ray tail 
and its evolution during long-term flux variations, ii) the variation ( 7 10 sec) of the best-fit 
spectral parameters for short- and long-term phases, and iii) the dependence of the spectral 
index and the electron temperature on the total flux and accretion rate. Data from the 
PCA and HEXTE detectors as well as BeppoSAX detectors have been used to constrain 
spectral fits, while ASM data provided long-term intensity state monitoring. Results of our 
long-term study of 4U 1820-30 are present, in detail, in the next sections and compared with 
our previous results for 4U 1728-34 and GX 3+1. 

We use the broadband energy spectra of BeppoSAX (Boella et al. 1997) and RXTE 
(Bradt, Rothschild, & Swank 1993) combined with the high-timing resolution of RXTE to 
study short and long term spectral and timing evolution of atoll sources. 
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3. Spectral Analysis 

Unlike the past analyses of the source spectral data discussed in the Introduction, in 
our study we make use of an unified model for both soft and hard states. In this way we 
have an opportunity to compare X-ray spectra of 4U 1820-30 in all states. 

In our spectral model, we use an assumption, that the accretion material passes through 
the accretion disk [for example, through the standard Shakura-Sunyaev disk (Shakura & Sunyaev 
1973)] and the transition layer (TL) (Titarchuk et al. 1998) where soft photons coming from 
the disk and NS surface are Comptonized off hot plasma (see also Fig. 2 in ST12). The Earth 
observer can also observe directly some fraction of these disk and NS seed photons. 

Thus, our input model is a sum of Comptonization component ( COMPTB ), which is 
the XSPEC Contributed model 3 , [see Farinelli et al. (2008), hereafter F08] and soft black- 
body and line ( Gaussian ) components. The parameters of the COMPTB component are 
the seed photon temperature kT s , the electron (plasma) temperature kT e: the energy index 
a (= r — 1 ) of the Comptonization spectrum, the illumination fraction of the Comptonized 
region (TL), / [/ = A/(l + A)] and the normalization of the seed photons illuminating the 
Comptonized region, N C omptb • We include a Gaussian component in the model charac- 
terized by the parameters Eu ne , <Ju ne Nu ne which are a centroid line energy, the line width 
and the normalization correspondingly. We also include a blackbody component and the in- 
terstellar absorption in our model characterized by the parameters: the normalization N BB , 
the color temperature T BB and a column density N H respectively. 

We fix the index of the seed photon spectrum a = 2 (or 7 = a + 1 = 3). Namely, we 
suggest that this seed photon spectrum is a blackbody-like. We neglect the bulk inflow effect 
with respect to the thermal Comptonization assuming that a bulk parameter <5 = 0. The 
parameter log(A) of the COMPTB component is fixed at 2 because the best-fit log(A) 1. 
Then / = A/(l + A) as the illumination fraction parameter is approximately 1, for any 
log(A) 1. We use a value of N H = 3.00 x 10 21 cm -2 estimated by Bloser et al. (2000) 
for 4U 1820-30. We find satisfactory fits using our model for both BeppoSAX. and RXTE 
observations of 4U 1820-30 for all available data sets. 


3 http: / /heasarc. gsfc.nasa.gov/docs/software/lheasoft/xanadu/xspec/models/comptb.html 


3.1. BeppoSAX data analysis 


Table 3 shows the data analysis results for the broad-band BeppoSAX spectra. On the 
top of Figure 1 we present an example of the BeppoSAX spectrum along with its best-fit 
using our model while in the bottom panel we demonstrate Ay (reduced y 2 =l.ll for 364 
d.o.f). The line emission is clearly centered around 6.7 keV. We find that the width of this 
line of 0.8 keV is quite large and it is much wider than the instrumental response which width 
is smaller than 0.02 keV 4 . This broad emission line at 6.7 keV can be a result of illumination 
of highly ionized iron by X-ray continuum. Piraino et al. (2000) suggest that this broad line 
originates either in an ionized innermost disk region or in a hot corona above the disk. A 
combination of two absorption edges related to ionized iron, instead of a Gaussian line, can 
also describe this part of the spectrum (D’Ai et al. 2006). The Laor relativistically smeared 
line or reflection models can be also used to describe this line feature (see Ng et al. 2010 
and Egron et al. 2011 respectively). However, Seihna & Titarchuk (2011) demonstrate that 
the model, wabs*(blackbody±COMPTB±Gaussian), which includes a a Gaussian iron line, 
can successfuly fit the data for extensive RXTE and BeppoSAX observations of 4U 1728-34. 

We can interpret this broad iron line detected in 4U 1820-30 in terms of reprocessing 
emission by a disk. In addition, we fit a smeared absorption edge in the 7—8 keV range using 
the smedge XSPEC model (see Ebisawa et al. 1994). The edge energy is 7.7±0.5 keV for 
the wabs * ( Bbody ± Comptb ± Gaussian) * smedge model which indicates to the presence of 
ionized material in the emission region. The smearing width is fixed at 10 keV. We include 
this edge component in the fits for all BeppoSAX data (see Table 3). 

We obtain that a = 1.03 ±0.04 (or T = a+ 1 = 2.03 ±0.04) for all analyzed BeppoSAX 
data, the seed photon temperature kT s of the COMPTB component is low variable and 
its value is around 1.3 keV whereas kTsB varies in the interval from 0.58 to 0.69 keV (see 
Tables 1, 3 for details). 


3.2. RXTE data analysis 

For all RXTE fits we fix the blackbody temperature kTsB =0.7 keV which is an upper 
limit of that in our analysis of the BeppoSAX data (see Table 3), because RXTE detectors 
cannot give us reliable spectra below 3 keV. In Table 4 we show the best-fit parameters of the 
RXTE spectra using our model. It is important to point out that for all RXTE observations 
of 4U 1820-30 the photon index T only slightly changes around 2 (T = 1.99 ±0.02) while the 


4 See ftp://heasarc.gsfc.nasa.gov/sax/cal/responses/98_ll 
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best-fit kT e varies in the 2.5—21 keV range. 

However the determination of the iron line profile using the RXTE data is a difficult 
problem caused by the low-energy resolution of PCA /RXTE detector. Moreover, the in- 
clusion of a smedge compoment in the spectral model for RXTE data does not improve fit 
quality any more. Therefore we apply our spectral model to RXTE data using a simple 
Gaussian as the line component without smedge modeling. The line width an ne does not 
vary much and it is always in the interval from 0.9 to 1.3 keV during all spectral transitions. 
Therefore we fix <Ju ne at 1.2 keV for all spectra during fitting procedure. The values of the 
best-fit seed photon temperatures, kT s = 1 — 1.3 keV are consistent with that obtained using 
the BeppoSAX data (see Table 3). 

In Figure 2 we show the representative examples of E * F(E) spectral diagrams of 
4U 1820-30 during soft ( left panel) and hard ( right panel ) state events. The best-fit RXTE 
spectra ( top panels) in the model wabs * ( Bbody + CompTB + Gaussian ) with Ay ( bottom 
panels) for the high- luminosity ( banana ) state [40017-01-11-00 observation, y^ ed =1.00 for 
78 d.o.f, left panel) and for the low- luminosity ( island ) state [94090-01-04-00 observation, 
Xred = 1-10 for 78 d.o.f, right panel). The model best-fit parameters are T = 1.99 ± 0.02, 
kT e = 2.94 ± 0.01 keV and E Gauss = 6.53 ± 0.06 keV for the soft state; T = 2.00 ± 0.04, 
kT e = 12.54 ± 0.09 keV and E Gauss = 6.35 ± 0.04 keV for the hard state (see Table 4 
for details). Violet, blue, red lines correspond to the Gauss, Blackbody and CompTB 
components, respectively. 

In Figure 3 to illustrate the spectral evolution 4U 1820-30 we show six representative 
EFe spectral diagrams for different electron temperatures of a Compton cloud [kT e =2.9 
keV (red), 3 keV (blue), 4 keV (green), 6 keV (violet), 10 keV (pink) and 12 keV (black)) 
applying the wabs * (Blackbody + COMPTB + Gaussian) model during island — banana 
state transitions. 

We show how the TL electron temperature kT e anti-correlates with the normalization 
IVcomptb (proportional to M) in Figure 4. For a comparison we add the points for 4U 1728-34 
and GX 3+1 [see Seihna & Titarchuk (2011) and Seihna & Titarchuk (2012) respectively]. 
The electron temperature kT e decreases and reaches a saturation about 3 keV when mass 
accretion rate increases [see an explanation of this effect in Farinelli & Titarchuk (2011), 
hereafter FT 11], 

Our spectral model applied to the spectral data of BeppoSAX and RXTE shows a robust 
performance throughout all data sets. Namely, a value of reduced xled = X 2 /Nd.o.f, where 
N d o f is a number of degree of freedom, is around 1.0 for most of observations, xled i s about 
1.5 for less than 3% of spectra with high counting statistics but xled never above a rejection 
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limit of 1.6. Note the high residuals of the poor fit spectra (2 among 234 spectra for which 
y 2 =1.55) occur in the iron line region. As was shown by BeppoSAX analysis a shape of the 
iron line is more complex than a simple Gaussian (see discussion in §3.1). Probably, the fits 
of this line indicate to a broad line, which shape and width can be a result of scattering of 
the line photons in the hot plasma (TL) along with iron smedge effect there. However, we 
cannot resolve this line complexity using the RXTE data. 

Thus using broad band BeppoSAX observations we can obtain the best-fit parameters 
of our spectral model whereas due to large-time coverage of 4U 1820-30 by RXTE we are 
capable to study the source spectral transitions in the 3—200 keV energy range. 


4. Overall pattern of X-ray properties 
4.1. Hardness-intensity diagram 

To study the properties of 4U 1820-30 during the spectral transitions when the luminos- 
ity changes we use hard color (10-50 keV/3-50 keV) (HC) versus the 3-60 keV flux measured 
in units of 10” 9 erg s -1 cm -2 [hardeness-intensity diagram (HID)]. 

In Figure 6 in TSF13 we demonstrate flux ratio HC versus the 3-60 keV flux using the 
RXTE data. As it appears from this Figure, 4U 1820-30 shows a “J” like diagonal shape in 
this diagram with upper short and lower elongated branches, which are joined at the lowest 
flux point. The spectral branches are indicated for the island state (IS), the lower banana 
state (LB) and the upper banana (UB) state. The direction of IS— >LB— ?-UB transition is 
shown by an arrow. The electron temperature kT e changes from 21 keV to about 3 keV 
along the arrow line direction (compare with Fig. 4). 

The identification of hardness- intensity diagram (HID) states are made using simulta- 
neous timing and spectral analysis, and we have revisited the previous similar RXTE data 
analysis made by Bloser et al. (2000) and Migliari et al. (2004). 

In particular, the hard color (HC) drops from 0.25 to 0.1 while the 3-60 keV flux is 
quasi constant when the source propagates from IS towards the LB. On the other hand HC 
rises from 0.12 to 0.35 with a simultaneous growth of the 3-60 keV flux when the source goes 
further from LB towards UB. 
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5. Evolution of X-ray spectral properties during spectral state transitions 

A number of X-ray spectral transitions of 4U 1820-30 with luminosity variations have 
been detected by RXTE during 1996 - 2009 (Rl - R7 sets). We investigate common 
spectral— timing signatures which can be found for these spectral transition events. The 
source reveals different behaviors during high-luminosity and low- luminosity events. 


5.1. Evolution of X-ray spectral properties during high-luminosity events 

In Figures 5—6 we show the results of our spectral analysis of the RXTE observations 
applying the wabs * ( blackbody + COMPTB + Gaussian ) model. On the top panels we 
present the RXTE/ ASM count rate and the model fluxes from 3 to 10 keV and from 10 to 
50 keV (see blue and crimson points respectively). The TL electron temperature kT e as a 
function of time is shown in third panel from above. The temperature kT e changes in the 
2.9—4 keV interval during the time period MJD 50490— MJD 50700 and only slightly varies 
around 3 keV during MJD 51200— MJD 51500 time interval. The COMPTB normalization 
N comptb and the normalization of low temperature blackbody component Nbb ( crimson 
and blue points respectively) are presented in the next-to-bottom panel of Figs. 5-6. One 
can clearly see that the COMPTB normalization Ncomptb correlates with variations of 
the ASM count rate and the 3-10 keV model flux. On the other hand the blackbody 
normalization Nbb only slightly varies. 


5.2. Evolution of X-ray spectral properties during the low-luminosity events 

Since the late April of 2009, 4U 1820-30 showed a prominent X-ray low-hard state at 
energies less than 10 keV as it was observed by the X-ray monitors on RXTE and Swift. We 
display the characteristics of the low-hard state obtained using the RXTE data in Figure. 7. 
Since 24 April 2009 (MJD 54945), the source was steadily brightening in the 15 - 50 keV 
band of the Swift/Burst Alert Telescope (BAT), with a daily average 0.032±0.002 cts/cm 2 /s 
(145 rnCrab) (Krimm et al. 2009). In contrast, the highest count rate detected by the 
Swift/ BAT was 0.14 cts/cm 2 /s. 

The ASM/ RXTE and PCA/ RXTE light curves showed that 4U 1820-30 were in an 
extended low state from MJD 54944 to MJD 54982. The ASM count rate took a sharp drop 
at MJD 54944 while the flux began rising in the BAT monitor. The ASM count rate was 
very low, approximately 6.0±0.5 cts/s during this low state period, with respect to an usual 
average count rate of ~20 cts/s. During the same time period the RXTE/PCA count rate 
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decreases from ~4000 cts/s to ~1000 cts/s. This kind of the long-time low state was not 
observed from 4U 1820-30 over the 10 - 15 year period. The typical low state duration varies 
from 1 to 2 weeks. 

We also establish that the X-ray spectra of this source over the low luminosity state 
(MJD 54955 - 54982) are quite stable in terms of the CompTB normalization value. But 
the electron temperature kT e of the Comptonised plasma increases from 3 keV up to 20 
keV during MJD 54950 - 54960 period and after that kT e gradually decreases again to 3 
keV when the luminosity rises at the end of the quasi-plateau (at MJD 55000). In Table 4 
we report the best-fit parameter values. During the IS-B transition period we find that the 
photon index T (or the spectral index a = T — 1) is almost constant, i.e. only slightly varies 
around 2 (see the combined Figure 10 in TSF13). 


5.3. Spectral state transitions in 4U 1820-30 

In 4U 1820-30 the hard-soft state transition is observed as kT e decreases from 15 - 20 
keV to 3 keV (see Fig. 7). The outburst hard-to-soft state transitions are seen when a 
supply of the soft photons flux N C omptb dramatically increases (see Fig. 4 and Fig. 6 in 
TSF13). In general, following FT11, we consider the spectral state transitions in terms of the 
kT e change. Thus the hard state is seen when the electron temperature reaches maximum, 
kT^ax , whereas the soft state is observed when the electron temperature, T e rmn reaches 
minimum. Note that kT e is well determined by the high energy cut-off of the spectrum, E cut 
and can be well established by the spectral fits to the data. This spectral state definition is 
based on kT e (or E cut ~ 2 kT e ). Unlike what occurs in the case of NS binaries, for a black 
hole case one can relate a spectral state change to the spectral (photon) index change [see 
Shaposhnikov & Titarchuk (2009), hereafter ST09]. 

We test the hypothesis of T appr ~ 2 using y; 2 -statistic criterion. We calculate the distri- 
bution of XredO^appr) = jj YliLi ( j versus of Y appr . and we find a sharp minimum of 
function Xredi^appr) at 1 when T appr = 1.99 ± 0.01 and Y appr = 1.99 ± 0.02 with confidence 
levels of 67% and 99% for 234 d.o.f. respectively (see the related Figure of XredO^appr) f° r 4U 
1728-34 in ST11). The index T is almost constant when kT e (see Fig. 8) and the COMPTB 
normalization, T39 /d\ Q changes (see below). FT11 based on the analysis of BeppoSAX data 
propose that T is about 2 for quite a few NS sources. FT 11 also define the spectral state 
using a value kT e and they demonstrate that T = 2 ± 0.2 (or a — 1 ± 0.2) when kT e varies 
in the 2.9—21 keV interval. 

It should be noted that not all NSs shows spectral state transitions but quite a few NSs 
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exhibit them. For instance, so called atoll-sources (such as 4U 1820-30), usually demonstrate 
IS - B transitions. Particularly during such transitions, it can be possible to differ NS from 
BH. Specifically, NSs and BHs show drastically different variations of spectral characteristics. 
NSs, as examples of 4U 1820-30 and 4U 1728-34, indicate to variabilities of kT e and mass 
accretion rate M along with a quasi-constant index T about 2. Meanwhile, BHs demonstrate 
a monotonic growth of V when M increases and succeeded by its final flattening (saturation) 
[see ST09]. 


6. Spectral— timing correlations during spectral state transitions 

We analyze the RXTE light curves applying the powspec task taken from FTOOLS 
5.1. We implement the timing analysis RXTE/PCA data which we perform in the 13-30 
keV range applying the event mode with time resolution of 1.2xl0“ 4 s. We make power 
density spectra (PDS) in 0.1—500 Hz frequency range with 0.001— second time resolution. 
The Poissonian statistics contribution was subtracted. We apply QDP/PLT package 5 for 
PDS modeling. 


6.1. Spectral and timing properties during low luminosity state transition 

In Figure 12 in TSF13 we present a generic behavior of X-ray timing— spectral char- 
acteristics for the low luminosity state at R7 (2009) transition events. We plot PDSs ( left 
column) along with the EF(E) — spectral diagram ( right column) for six moments at MJD = 
54947.6/54956.5, 54958.6/55002.6 and 54997.7/54960.36, covering different transition phases. 
At the bottom we demonstrate PDSs for the 15—30 keV energy range ( left column) and plot 
along with the E * F(E) — spectral diagram ( right column) for A— C time events (see upper 
panel). All points [(events) A red (ID 94090-01-01-00), A blue (ID 94090-01-02-03), B red 
(ID 94090-01-02-02), C red (ID 94090-01-04-03)], except B (blue) and C (red), are related 
to IS (broadband noise, no VLFN). 

PDSs denoted by B (blue, ID 94090-02-02-00) and C (red, ID 94090-0103-00) exhibit the 
island— lower banana state transition. For the blue PDS the VLFN (very low frequency noise) 
appears as a band-limited noise component which transforms into QPO (broad Lorentzian 
with u h centroid frequency at 7 - 10 Hz). We present power density spectra (PDSs) (panels 
Al, Bl, Cl) along with the corresponding E * F (E)— spectral diagrams (panels A2, B2, C2). 


5 http: / /heasarc. gsfc.nasa.gov/ftools/others/qdp/qdp.html 
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The related PDSs and energy spectral data are shown by blue and red points respectively. 
On the left panel we also show the electron temperature kT e associated with a given PDS. 


6.2. Spectral and timing properties during high luminosity state transition 

To compare with Figure 12 in TSF13, we show an evolution of spectral— timing charac- 
teristics during the high luminosity state in Figure. 9. Here on the top panel we display the 
ASM light curve during high luminosity interval at R3 (1999) transition events. Red/blue 
points A, B, and C are related to the moments at MJD = 51283.6/51300, 51313.7/51330.5 
and 51389.4/51396.26, covering different transition phases. On the bottom left and right 
panels we present PDSs for 15-30 keV energy range and the E * F(E) — spectral diagram cor- 
respondingly, for A (red, top), B (blue, middle) and C (blue, bottom) points of X-ray light 
curve (see the upper panel). All points are related to the banana state with relatively strong 
broadband noise and VLFN with QPOs which is at ~ 6 - 7 Hz for C moment (red). The 
PDSs in panel C (blue) and (red) illustrate the island— lower banana state transition. Here 
the VLFN appears, band-limited noise component transforms into QPO [a broad Lorentzian 
with Uh centroid frequency at 7—10 Hz, C red (40017-01-12-00)]. The power spectra (panels 
Al, Bl, Cl) correspond to E * F(E) diagrams (panels A2, B2, C2). The corresponding 
energy spectra of 4U 1820-30 are related to the electron temperature of 3 keV. 

In Figure 10 we illustrate a typical power spectrum of 4U 1820-30 for different X-ray 
spectral states (shown on the right panel). The electron temperature values of corresponding 
energy spectra are indicated at the right vertical axis. The power spectra in the extreme 
island state (EIS), island state (IS, multiplied by factor 10 -2 for clarity), lower left ba- 
nana state (LLB, xl0~ 4 ), lower banana state (LB, xlO -6 ) and upper banana state (UB, 
xl0~ 8 ) are presented from the top to the bottom. The histograms show the best fits to 
the power spectra, which consist of three components: VLFN the peaked noise component, 
low-frequency QPOs (zy and Uh) and high frequency feature VhHz (see van Straaten, van der 
Klis & Mendes 2003 for details of terminology). 


6.3. Comparison of spectral and timing characteristics of atoll sources 

4U 1820-30, GX 3+1 and 4U 1728-34 

In this Paper, we also study the correlations of X-ray spectral— timing characteristics 
and M in a number of atolls during their spectral transitions searching for similarities and 
differences between atoll sources. In this way we can present a comparative analysis for 
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three atoll sources: 4U 1820-30, GX 3+1 and 4U 1728-34 applying the same spectral model 
which consists of low temperature Blackbody , Comptonized continuum and Gaussian line 
component. 


6.3.1. Constancy of the photon index 

We demonstrate that atolls 4U 1820-30, GX 3+1, 4U 1728-34 show a similar pattern of 
the photon index Y vs M (or Ncomptb)- Namely, the photon index Y only slightly varies 
around 2 (see Fig. 11). Following to FT11, ST11 and ST12, we can suggest that the cooling 
flow of soft disk photons is much less less than the energy release in the transition layer (TL) 
for each of these three sources. 


6.3.2. The difference and similarity of kT e ranges in fU 1820-30, GX 3+1 and 

4U 1728-34 

One can see from Figure 12, that the ranges of kT e for an individual state evolution of 
these three sources are different. The electron temperature kT e changes in 4U 1728-34 from 
3 to 15 keV whereas kT e varies within much narrow range of kT e around 3 keV in GX 3+1. 
In turn, source 4U 1820-30 demonstrates a wider interval of kT e in which kT e varies from 2.9 
keV to 21 keV similar to some extent to the temperature change in 4U 1728-34. Note, that 
in a low temperature regime 4U 1820-30 and GX 3+1 are similar in terms of normalization 
-Ncomptb =(4 - 15) xL 3 g/D'f 0 , or mass accretion rate, (see Fig. 4) and Comptonized fraction 
/ = 0.2 - 0.8 (see Fig. 12). While 4U 1820-30 and 4U 1728-34 are similar for a range of 
the normalization Womptb =(2 - 4fxL 39 / Df 0 (see Fig. 4) and / = 0.5 — 0.8 (see Fig. 12) for 
high temperature regime. However, in contrast to 4U 1728-34, the source 4U 1820-30 has an 
additional branch of intermediate temperatures (8 - 12 keV) when the Comptonized fraction 
is relatively low, / < 0.5 ( see Fig. 12). Note that all objects have a common temperature 
interval 3-4 keV when Womptb (or mass accretion rate) is relatively high. 

Thus according to FT11, ST11, ST12 and the present study, the electron temperature 
kT e for atolls and Z— sources varies in the 2.5—25 keV range. Specifically, the change of 
kT e around 3 keV is similar for all three source 4U 1820-30, GX 3+1 and 4U 1728-34. The 
minimum value of kT e at 2.5 keV occurs at the peak luminosity for 4U 1728-34 (see ST11), 
during a local rise of luminosity for GX 3+1 (ST12) and at high luminosity phases for 
4U 1820-30. 

For all of these three objects, the values of color seed photon temperatures kT s = 1.1 — 1.7 
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keV and blackbody temperatures kT BB ~0.6 keV are comparable (see Table 5). Contrary, the 
variability extent of kT e is not similar. The reason for that difference of electron temperature 
ranges is evident. Sources 4U 1820-30 and 4U 1728-34 show a complete cycle of state 
evolution: island— lower banana (LB)- upper banana (UB) stages for 4U 1820-30 and for 
4U 1728-34 it is extreme island state (EIS )— upper banana (UB) state (see Di Salvo et al. 
2001; ST11). But GX 3+1 demonstrates a short evolution behavior on the CCD from LB 
to UB. This evolution picture is also clear from Figure 12 which shows that the track of 
GX 3+1 is only a part of the full track [see definition of a state sequence and the standard 
atoll- Z scheme in Hasinger & van der Klis (1989)]. 

Note, that 4U 1820-30 shows almost the same kT e range as that in 4U 1728-34 and 
almost similar timing evolution. But clear differences between these atolls one can see from 
Figure 17 in TSF13 where we show spectral hardness (10 - 50 keV/3 - 50 keV) vs flux in 
3-60 keV range. In fact, 4U 1728-34 ( blue points) is fainter and harder and demonstrates 
much wider spectral hardness range than that in 4U 1820-30 (red points). 


6.3.3. Comparison of spectral evolution as a function of the luminosity for fU 1820-30, 

GX 3+1 and fU 1728-3 f 

Now we present a comparison of X-ray spectrum evolution for sources 4U 1820-30, 
GX 3+1 and 4U 1728-34 based on luminosity value which is presumably proportional to 
Comptb normalization and, consequently, to mass accretion rate taking into account that 
their distances to the Earth are similar (see Table 5). For 4U 1820-30 the distance range is 
within of 5.8 - 8 kpc (Shaposhnikov & Titarchuk, 2004), whereas for 4U 1728-34 and GX 3+1 
the distances are estimated as 4.5 kpc and 4.2— 6.4 kpc respectively [see van Paradijs (1978) 
and Kuulkers & van der Klis (2000)]. 

We show the CompTB normalization (related to the soft photon luminosity value) for 
these sources as a function of kT e in Figure 4. 4U 1820-30 subtends a wider interval in 
CompTB normalization than that for 4U 1728-34. Note, that in the high luminosity state 
(or Womptb) 4U 1820-30 is similar to GX 3+1: Womptb = (4 - 15)xL 3Q / Df 0 , Comptonized 
fraction / =0.2 - 0.8 (see Fig. 12) and the electron temperatures kT e are low variable around 
3 keV. While in the low luminosity state 4U 1820-30 is closer to source 4U 1728-34: !V C omptb = 
(2 - 5 )xL 39 /Df 0 (see Fig.4), / =0.5 - 0.8 and kT e changes from 5 to 20 keV (see Fig. 12). 
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6.3.4- The difference and similarity of time scales of state evolution for 4 U 1820-30, 

GX 3+1 and 4U 1728-34 

We should point out that all these three atoll source show the transitions between 
low luminosity and high luminosity states over different time scales. Specifically, the time 
scales of X-ray flux variability for 4U 1728-34, 4U 1820-30 and GX 3+1 probably dictated 
by variability of mass accretion rate, are ~10 d, 100 d and 1000 d, respectively. However, 
these sources demonstrate LB - UB transition and make it in the narrow interval of low 
temperature kT e (around 3 keV) and during the same short time interval (hours - day). We 
remind the reader that the comparison between these three sources is facilitated by the fact 
that they show almost the same kTsB and kT s temperature values and they are located at 
approximately the same distance. The only difference of spectral evolution of these objects 
is related to different ranges of the electron temperature of the Comptonized component. 


6.3.5. Correlation of illumination parameter f versus electron temperature kT e and its 
relation with different stages in the color-color diagram 

Using Table 5 one can see that the ranges of the best-fit illumination fraction / are 
0.2 — 1.0, 0.2 — 0.9 and 0.5 — 1 for 4U 1820-30, GX 3+1 and 4U 1728-34 correspondingly. 
These values of / indicate to different geometry of the transition layer (TL) and thus to 
different illumination for these X-ray sources. 

In Figure 12 we demonstrate that the electron temperature kT e directly correlates with 
a sequence of CCD states, EIS-IS-LLB-LB-UB [see Hasinger & van der Klis (1989) for this 
CCD classification]. Note that ST12 reveal a relation between spectral states, kT e and / for 
atoll sources GX 3+1 and 4U 1728-34. We show these kT e — f relations for these two atolls 
in Figure 12. The direction in which the inferred M increases is indicated by arrows. 

Now we present three different tracks on the kT e — / diagram for three source 4U 1820- 
30, GX 3+1 and 4U 1728-34 and show how these tracks are related to the standard CCD 
sequence (see Fig. 12). The track of 4U 1820-30 consists of three segments (branches) related 
to kT e : high (12 - 21 keV), intermediate (7 — 12 keV) and low temperature (2.9 - 6 keV) 
ones, wherein each segment has a negative correlation of kT e and /. In turn, GX 3+1 
demonstrates the only, so called, low temperature branch track. Namely, when the fraction 
/ increases, kT e decreases from ~ 4.5 keV to ~ 2.3 keV. While for 4U 1728-34 we see a 
more complicated pattern, but in contrast to 4U 1820-30 and GX 3+1, has a segment with 
positive correlation of kT e vs / from 4 to 12 keV. Specifically, at the high temperature state 
(EIS), / only slightly changes from 0.9 to 1 when kT e decreases. As kT e further drops from 


12 keV to 4 keV, / also drops from 0.9 to 0.5. Finally, / goes up from 0.5 to 0.8 when the 
source enters to the low-temperature state (LB-UB). As a result, we demonstrate that the 
CCD state evolution can be also seen using the kT e — f correlation. 


6.4. Comparison of spectral hardness diagrams for atolls 4U 1820-30, GX 3+1 

and 4U 1728-34 

We use the plot HC (10-50 keV/3-50 keV) versus the 3-60 keV flux in form of HIDs 
for three sources: 4U 1820-30 (red), GX 3+1 (green) and 4U 1728-34 (blue) (see Fig. 17 in 
TSF13) to compare transition properties of these atolls in terms of their flux (or luminosity). 
In fact, 4U 1820-30 shows a “J” like inclined (or diagonal) shape in HID with upper short 
and lower elongated branches (see Fig. 6 in TSF13). The short branch is close to the low 
luminosity state, whereas the elongated branch covers the wide luminosity range. 

Our comparative analysis of HID track branches for 4U 1820-30, GX 3+1 and 4U 1728-34 
indicates to similar physical properties of these objects. The spectral and timing character- 
istics are very similar along corresponding segments. Specifically, the short branch of a “J” 
like track of 4U 1820-30 is adjacent to low luminosity area of 4U 1728-34 (blue points) and 
it is related to a high electron temperature regime of 4U 1820-30, as in 4U 1728-34. In turn, 
the elongated branch of 4U 1820-30 is closer to GX 3+1 (green points) branch area and it 
is associated with low electron temperatures kT e (3-4 keV) and softer spectra, which are 
also seen in GX 3+1. 

Note that among considered atolls superbursts are only observed in GX 3+1 and 
4U 1820-30 during elongated branch. Furthermore, superbursts are detected at low— soft 
states, i.e. during low luminosity interval of light curve when electron temperature kT e is 
about 4 keV. Thus 4U 1820-30 shows a property similar to GX 3+1 and it is situated at 
an intermediate position between 4U 1728-34 and GX 3+1 in terms its luminosity. This 
observational fact can be related to the same intermediate rate of mass transfer in these two 
sources, 4U 1820-30 and GX 3+1 [see also a review by van der Klis (1994)]. The comparison 
of HIDs allows to diagnose physical properties of different objects with adjacent HID tracks. 
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7. Discussion 

7.1. Stability of photon index is a signature of NS source 

Thus we demonstrate that the photon spectral index only slightly varies around 2 using 
numerous observations of NS sources 4U 1820-30, GX 3+1 (ST12) and 4U 1728-34 (ST11) 
by BeppoSAX and RXTE. In Figures 8, 11 and 13 ( left column) we show T as a function 
of the spectral model parameters: kT e (in keV), Ncomptb— normalization, and illumination 
fraction /. These results for NS 4U 1820-30 have been obtained when we apply our thermal 
Comptonization model to BeppoSAX. and extensive RXTE observations. FT 11 and ST 11 
also investigate and find the photon (energy) index stability in other observations of NS bi- 
naries. We explain this index stability using the Comptonization model. Namely, the photon 
(energy) index is almost constant when the soft photon flux illuminated the transition layer 
(TL) is much less than the gravitational energy release in TL (see e.g. ST12). This model 
of the index stability can probably resolve the index stability effect now clearly established 
in these three NS sources using extensive BeppoSAX and RXTE observations. 


7.2. On the hard tail origin in atoll source 4U 1820-30 

The radio emission detected from 4U 1820-30 (Migliari et al. 2004) suggests the pres- 
ence of a jet, which may also generate an extended power-law X-ray emmision. In this 
case, the power law can be a result of the inverse Compton effect on nonthermal electrons 
of the jet. Note, X-ray nonthermal power-law tails are also observed in soft states of BHs 
see for example a review by ST09; see also McConnell et al. (2002) and Wardzi’nski et al. 
(2002) on the detection of the extended hard tails in the hard states of BHs, Cyg X-l and 
GX 339-4 respectively, and NS Z-sources [see Di Salvo et al. (2000); Farinelli et al. (2005); 
D’Amico et al. (2001); Asai et al. (1994)]. However these extended hard tails are also 
found in atolls [see e.g. Piraino et al (1999)]. 

Additive models that have been applied to fit the spectra of 4U 1820-30, need to use 
an additional power-law component (pure one or as a component of CompPS) to describe 
a hard spectral tail above 80 keV [see e.g. Tarana et al. (2006)]. However such an ap- 
proach invokes an unknown non-thermal origin of hard tail emission. On the other hand our 
suggestion allows us to explain X-ray spectra of 4U1820-30 in all spectral states using the 
same model without a specific composition of the model components at different states. In 
fact, in our model (see §5.2 and §6.1) we describe the hard tail emission using the thermal 
Comptonization component in which the TL electron temperature kT e increases up to 20 
keV and the illumination factor / decreases as the source goes to the hard state (see Figs. 
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3, 12). 


8. Conclusions 

We analyze the X-ray spectral and timing characteristics of 4U 1820-30 observed during 
the hard-soft state transitions. We find a number of spectral transitions in 4U 1820-30 using 
BeppoSAX and RXTE data. 

For our investigation we take an advantage the BeppoSAX broad spectral extension over 
the 0.3—200 keV range and abundant RXTE observations taken in the 3—200 keV energy 
coverage. 

We demonstrate that the X-ray broad-band spectra can be successfully fit by com- 
position of the Blackbody , Comptonization ( Comptb ) and Gaussian— line components for all 
spectral states. Also we show an observable relation of the photon index T and the normaliza- 
tion of the Comptonized component, Comptb which is proportional to M. We demonstrate 
the stability of the photon index T ~ 2 when the source goes from the hard state to soft state, 
in other words when the electron temperature of Comptonized region (TL), kT e decreases 
from 21 to 3 keV (see Fig. 8). 

We also show that T only slightly varies with the Comptb normalization (oc M). Note, 
this stability of the index in NS sources has been recently suggested for a number of other 
NSs, Sco X-l, Cyg X-2, GX 17+2, GX 3+1, GX 340+0, GX 349+2, X 1658-298, IE 1724- 
3045, GS 1826-238, which were analyzed using BeppoSAX data (see details in FT11, ST11, 
ST 12). The use of the disk seed photon normalization, (Comptb), which is proportional to 
M, is fundamental in order to find the stability of F during the hard-soft state transition. 
We do find the stability (constancy) of the photon index of Comptonized component versus 
both the Comptb normalization and the electron temperature kT e about 2 for all spectral 
states. In our analysis of NS sources (see FT11, ST11, ST12, Seifina et al. 2013 and this 
Paper) we do not find any particular case in which the photon index T changes beyond the 
limits 2 ± 0.1. Thus this index stability can be taken as an intrinsic property of neutron 
star (NS) binaries (as a NS signature ), which is drastically different from that in black 
hole binaries [e.g., GX 339-4, GRO J1655-40, XTE J1650-500, XTE J1550-564, 4U 1543-47, 
XTE J1859+226, H 1743-322, (ST09), GRS 1915+105 (TS09), SS 433 (ST10)], where T 
monotonically rises during the hard— soft state transition and it follows by its saturation at 
high M— values (see Fig. 13). In Figure 13 we show the T — M correlation for a number of 
BHs (right column) and that T is almost independent of M in NSs (left column). Indices 
T in BHs show clear correlation with M , or with the normalization Lsg/Df 0 [where L39 is a 
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flux of soft (seed) photons]. The Y — M correlation is followed by T— saturation when mass 
accretion rate M exceeds the Eddington limit. The behavior of T vs M for a considered 
sample of NSs (4U 1820-30, 4U 1728-34 and GX 3+1) is drastically different from that for 
given examples of BHs. 

A relatively wide interval of the illumination fraction / = 0.2 — 1 which we obtain in the 
framework of our model, point to variable soft (disk) photon illumination of the transition 
layer in 4U 1820-30. Using BeppoSAX data we also find two types of blackbody photons. 
One type is characterized by color temperature of 0.7 keV, which is typical for the disk 
photons and another one is related to 1.3 keV, which can be associated with NS surface 
temperature. 

We detect an evolution of 6 - 20 Hz QPOs and noise components during the island - 
banana state evolution (LLB-UB) (see Fig. 9). 

Our observational results establishing the constancy of the photon index Y in fU 1820- 
30 confirm the theoretical arguments of FT11 and ST11 that the TL energy release Q cor 
dominates the soft photon flux illuminating the transition layer (TL) which comes from the 
accretion disk, Qdisk ■ We argue that the found stability of Y is an intrinsic NS signature as 
in BH binaries Y monotonically increases with M followed by its saturation at high values 
of M (see ST09). 

LT acknowledges discussion with Chris Shrader and his thorough editing of the manuscript 
and also we appreciate comments of the referee which substantially improve a quality of the 
presented material. 
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Table 1. The list of BeppoSAX observations of 4U 1820-30 used in analysis. 


Obs. ID 

Start time (UT) 

End time (UT) 

MJD interval 

20105004 

1997 Oct. 2 06:45:23 

1997 Oct. 2 19:09:36 

50723.2-50723.8 

20537004 

1998 Apr. 17 04:31:23 

1998 Apr. 18 02:55:28 

50920. 1-50921. 1 1 

20537005 

1998 Sept. 23 12:44:56 

1998 Sept. 24 15:30:05 

51079. 5-51080. 6 1 


Reference. (1) Kaaret et al. (1999) 


Table 2. The list of groups of RXTE observations of 4U 1820-30 


Number of set 

Dates, MJD 

RXTE Proposal ID 

Dates UT 

Ref. 

R1 

50204-50207 

10074 

May 1 - 4, 1996 

1 


50151, 50235 

10076 

March 9, June 1, 1996 

1 


50371-50386 

10075 

Oct. 15 - 30, 1996 

1, 2, 3, 4 

R2 

50488-50701 

20075 

Feb. 9 - Sept. 10, 1997 

1, 3, 4, 5, 6, 7 

R3 

50920-51467 

30053, 30057 

Apr. 4, 1998 - Oct. 16, 1999 

3 


51206-51956 

40017, 40019 

Jan. 28, 1999 — June 4, 2003 


R4 

51996-51999 

60030 

March 28 - 31, 2001 



52355, 52429 

40017 

March 22, 2002; June 6, 2003 



52439 

70031 

June 14, 2002 13:21:52 - 14:14:40 



52482-52808 

70030 

July 23, 2002 - June 18, 2003 

8 


52894 

80105 

Sept. 12, 2003 05:41:04 - 09:36:32 

8 

R5 

53258-53591 

90027 

Sept. 10, 2004 - Aug. 9, 2005 



53692-53693 

91435 

Nov. 18 - 19, 2005 


R6 

53959-54028 

91151 

Aug. 12 and Oct. 9, 2006 



53921-54306 

92030 

July 5, 2006 - July 25, 2007 



54126-54129 

70030 

Jan. 26 - 29, 2007 


R7 

54947-55002 

94090 

Apr. 26 - June 20, 2009 

9 


54947-55116 

92030 

Sept. 8 - Oct. 12, 2009 



References: (1) Chou & Grindlay (2001); (2) Smale et al. (1997); (3) Kaaret et al. (1999); (4) Zhang et al. (1998); (5) 
Bloser et al. (2000); (6) Kusmierek et al. 2011; (7) Shaposhnikov & Titarchuk (2004); (8) Migliari et al. (2004); (9) Krimm 

et al. (2009) 


Table 3. Best-fit parameters of spectral analysis of BeppoSXX observations of 4U 1820-30 in 0.3-200 keV energy 

range'. Parameter errors correspond to 90% confidence level. 
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Table 4. Best-fit parameters of spectral analysis of PCA&HEXTE / RXTE observations of 4U 1820-30 in 3-200 keV 

energy ranged Parameter errors correspond to 90% confidence level. 
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Table 5. Comparison of parameters of atoll sources 4U 1820-30, GX 3+1 1 and 4U 1728-34 : 
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Fig. 1. — Top : the best-fit spectrum of 4U 1820-30 during banana branch events in E* F(E) 
units using BeppoSAX. observation 20105004 carried out on 2 October 1997. The data are 
presented by crosses and the best-fit spectral model wabs*(blackbody+Comptb+Gaussian) by 
green line. The model components are shown by bine, red and crimson lines for blackbdody, 
Comptb and Gaussian components respectively. Bottom panel. Ay (reduced y 2 =l.ll for 
364 d.o.f). The best-fit model parameters are T = 2.00 ± 0.04, kT e = 3.25 ± 0.02 keV, 
Eune = 6.7 ± 0.1 keV (see more details in Table 3). 
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Fig. 2. — Examples of typical E*F(E) spectral diagram of 4U 1820-30 during soft ( left panel ) 
and hard ( right panel ) state events. The best-fit RXTE spectra ( top panels) using the model 
wabs * ( Bbody + CompTB + Gaussian) with Ay ( bottom panels) for the high- luminosity 
(banana) state [40017-01-11-00 observation, y^ ed =1.00 for 78 d.o.f, left panel) and for the low- 
luminosity (island) state [94090-01-04-00 observation, y^ ed =1.10 for 78 d.o.f, right panel) . 
The best-fit model parameters are T=1.99±0.02, kT e = 2.94±0.01 keV and E Gauss =6.53±0.06 
keV (for the soft state) and T=2.00±0.04, kT e = 12.54±0.09 keV and E , Gauss =6.35±0.04 keV 
(for the hard state) (see more details in Table 4). Blue, red and violet lines stand for Bbody, 
CompTB and Gauss components respectively. 
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Fig. 3. — Six representative EFe spectral diagrams which are related to different electron 
temperatures of TL [kT e =2.9 keV (red), 3 keV (blue), 4 keV (green), 6 keV (voilet), 
10 keV (pink) and 12 keV (black)] using the model wabs * (Blackbody + COMPTB + 
Gaussian) for island — banana state transitions of 4U 1820-30. The data are taken from 
RXTE observations 30057-01-04-01 (red), 70030-03-07-020 (blue), 70030-03-05-02 (green), 
70030-03-05-01 (violet), 40017-01-24-00 (pink) and 94090-01-04-00 (black). 
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Fig. 4. — Comptb normalization measured in units of L 8 ^ 1 /D\ q versus the electron temper- 
ature kT e (in keV) obtained using the best-fit spectral model wabs * ( blackbody + Comptb + 
Gaussian ) for atoll sources 4U 1820-30 (red), GX 3+1 (green, taken from ST12) and 
4U 1728-34 (blue, taken from ST11) for RXTE data. 
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Fig. 5. — From Top to Bottom: Evolutions of the RXTEf ASM count rate, the model 
flux in 3-10 keV and 10-50 keV energy ranges ( blue and crimson points respectively), the 
electron temperature kT e in keV, and Comptb and blackbody normalizations ( crimson and 
blue respectively) during 1996 - 1997 transition set (R1 - R2 ). The rising phases of the local 
(mild) transitions are marked with blue vertical strips. 
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Fig. 6. — Similar to that presented in Fig. 5 but for the RXTE 1999 transition set R3. 
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Fig. 7. — Similar to that presented in Fig. 5 but for the RXTE 2009 transition set R7. 
The quasi-plateau phases of the low luminosity state of 4U 1820-30 are marked using orange 
vertical strip. 
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Fig. 8. — The photon index T plotted versus the electron temperature kT e (in keV) in the 
frame of our spectral model wabs * (blackbody + C omptb+ Gaussian) during transition events 
(see Tables 3, 4). Blue and red points correspond to BeppoSAX and RXTE observations. 
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Fig. 9. — Top-, evolution of RXTEf ASM count rate during the high luminosity 

state at R3 (1999) transition events. Red/blue points A, B, and C mark moments at 
MJD = 51283.6/51300, 51313.7/51330.5 and 51389.4/51396.26 covering different transi- 
tion phases. Bottom : PDSs for 15-30 keV band ( left column) are plotted along with the 
E * F(E )— diagram ( right column) for A (red, top), B (blue, middle) and C (blue, bottom) 
points of X-ray light curve. All points are related to the banana state [strong broadband 
noise, VLFN and QPOs at Vi ~6 - 7 Hz (C red)}. The E * F(E )— diagrams (panels A2, 
B2, C2) are related to the corresponding power spectra (panels Al, Bl, Cl). The data are 
shown by black points. The electron temperature kT e of the corresponding energy spectra 
of 4U 1820-30 is about 3 keV. 
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Fig. 10. — PDSs of 4U 1820-30 related to its X-ray spectral states. kT e values (in keV) corre- 
sponding to the energy spectra are indicated at the right vertical axis. PDSs in the extreme 
island state (EIS), island state (IS, multiplied by factor 10“ 2 for clarity), lower left banana 
state (LLB, xlO” 4 ), lower banana state (LB, xlO” 6 ) and upper banana state (UB, xlO” 8 ) 
are presented from top to bottom. The histograms consist of three components: VLFN ( very 
low frequency noise in banana states), the peaked noise component, low-frequency QPOs are 
fit by Lorentzians (ui, zy,,) and high frequency QPO ( VhHz )■ 
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Fig. 11. — From Top to Bottom: Plots P vs Comptb Normalization ( left column ) and vs 
Comptonized fraction / ( right column) for GX 3+1 (top), 4U 1820-30 ( middle ) and 4U 1728- 
34 ( bottom ) obtained using the wabs * ( blackbody + Comptb + Gaussian ) model. On the top 
panels crimson and blue points are for GX 3+1 taken from ST12 and on the middle panels 
for 4U 1820-30 red and blue points correspond to RXTE and BeppoSAX data respectively 
(current study). On the bottom panels blue and green points correspond to RXTE and 
BeppoSAX data respectively for 4U 1728-34 (data taken from ST11). 


50 



EIS 

IS 

LLB 

LB 

UB 


Fig. 12. — kT e (in keV) plotted versus illumination fraction / for 4U 1820-30, GX 3+1 (taken 
from ST12) and 4U 1728-34 (taken from ST11) during spectral state transitions obtained 
using the wabs * ( blackbody + Comptb + Gaussian ) model. Red, green and blue points 
correspond to RXTE observations of 4U 1820-30, GX 3+1 and 4U 1728-34 respectively. The 
bended arrows are related to an increase of mass accretion rate. On the right-hand side of 
the Figure we show a sequence of CCD states (EIS - extreme island state, IS - island state, 
LLB - lower left banana state, LB - lower banana state and LIB - upper banana state) which 
are listed according to the standard atoll— Z scheme (Hasinger & van der Klis 1989). One 
can see that kT e is directly related to the sequence of CCD stages. 
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Fig. 13. — Examples of diagrams of the photon index T versus the COMPTB normalization 
(proportional to mass accretion rate) for BH sources [right column, GRS 1915+105 (TS09), 
GRO J1655-40 (ST08) and GX 339-4 (ST08)] along with atoll NS sources [left column, 
4U 1820-30, 4U 1728-34 (ST11) and GX 3+1 (ST12)]. For all plots the RXTE data were 
used along with BeppoSAX data (indicated by blue points on the left column). One can 
see a noticeable change of T followed by the saturation plateau for BHs as for NSs the index 
only slightly varies about 2. The level for T = 2 is indicated by blue dashed line. 


